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Abstract 

The Berezin-Klauder-Toeplitz ( "anti-Wick" ) quantization or "non-commutative read- 
ing" of the complex plane, viewed as the phase space of a particle moving on the line, is 
derived from the resolution of the unity provided by the standard (or gaussian) coher- 
ent states. The construction properties of these states and their attractive properties 
are essentially based on the energy spectrum of the harmonic oscillator, that is on 
the natural numbers. This work is an attempt for following the same path by con- 
sidering sequences of non-negative numbers which are not "too far" from the natural 
numbers. In particular, we examine the consequences of such perturbations on the 
non-commutative reading of the complex plane in terms of its probabilistic, functional, 
and localization aspects. 

1 Introduction 

Ever since the initial discovery of coherent states by Shrodinger in 1926 and their redis- 
covery in the early sixties, independently by Glauber and Sudarshan (in the context of 
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paris7.fr, regomont@cbpf.br, ligia@cbpf.br, ken.sekimoto@espci.fr 
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quantum optics) and by Klauder (following a more thorough examination of the relation- 
ship between classical and quantum mechanics), many striking features and applications of 
these particular superpositions of harmonic oscillator eigenstates have been unveiled and 
explored. Certainly, one of the most intriguing and deep properties of these states is the 
Bayesian duality Ali, Gazeau, Heller (2008) | that they encode between the discrete Poisson 



probability distribution, n \— > e~^ 2 \z\ 2 /n\, of obtaining n quantum excitations ("photons" 
or "quanta" ) in a measurement through some counting device, and the continuous Gamma 
probability distribution measure \z\ 2 i— ► e - ' 2 ' |z| 2 /n! on the classical phase space, with z 
being a parameter defining the coherent states. For this latter distribution, \z\ 2 is itself a 
random variable, denoting the average number of photons, given that n photons have been 
counted. Such a duality underlies the construction of all types of coherent state families, 
provided they satisfy a resolution of the unity condition. It turns out that this condition 
is equivalent to setting up a "positive operator valued measure" (POVM) on the phase 
space. Such a measure, in turn, leads in a natural way to a quantization procedure for 
the classical phase space, namely the complex plane, a procedure which we also look upon 
as being a "non-commutative reading" of the complex plane. (The non-commutativity is 
understood in the sense that we associate to each phase space point z the one dimensional 
projection operator P z , projecting onto to the subspace generated by the coherent state 
vector, and then for z 7^ z' , P^Py 7^ P a /P a ). This "Berezin-Klauder-Toplitz" quanti- 
zation turns out, in this case, to be equivalent to the canonical quantization procedure. 
Clearly, this non-commutative reading of the complex plane is intrinsically based on the 
nonnegative integers (appearing in the n\ term). In this paper we will follow a similar 
path by considering sequences of non-negative numbers which are not "too far" from the 
natural numbers. The resulting quantizations will then be looked upon as perturbations of 
the non-commutative reading given by the standard coherent states. In particular, we will 
examine the consequences of such perturbations on the non-commutative reading of the 
complex plane as they relate to their probabilistic, functional, and localization properties. 

More precisely, let us consider an infinite, strictly increasing sequence of nonnegative 
real numbers {x n }neN such that xo = and there exist r, L € R such that < r < 
x n+ \ — x n < L < 00 for any n 6 N ("Delone sequence"). To this sequence of numbers 
correspond the sequence of "factorials" x n \ = x\x<i . . . x n with xq! = 1) the "exponential" 



m = E 

the sequence of "moment" integrals, 



00 f 



n=0 



r I ' 



-00 

X n-H"n = I J\f(t) ^ ' 







and the "renormalized" sequence 



x n := x n , n > 1 , x = . 
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In the present study we consider various examples of such sequences, for which it is 
proved, analytically or numerically, that the ratios — — ► 1. 

The "non-commutative reading" of the complex plane is understood in the following 
sense. Standard coherent states are vectors labelled by z G C, v z = e - ^ I 2 X] ne M z n /y/riJ.e n , 
in some separable Hilbert space with orthonormal basis {e n , n G N}. One of their most 
striking features is that they solve the unity: 

I = / (V z , -)Vz ■ 

Jc n 

This property allows us to transform ordinary functions f(z), z = (l/y/2)(q + ip), into 
f <Pz 

operators / h- > At = / f(z)(v z , -)v z ("coherent state quantization") and yields a non- 

Jc T 

commutative version of the plane illustrated by = il, where Q = A q and P = A p . 

Moreover, these states minimize that non-commutativity in the sense that the uncertainty 
relations (or Heisenberg inequalities) for the quantized coordinates Q and P are saturated: 



AQAP = - , AQ = y/(v z ,Q 2 v z )-(v z ,Qv z )*. 

Such an equality produces a "fuzzy" vision of the plane with elementary "minimal" square 
lattice cells of side length 1 / y/2. 

Deformed coherent states are constructed using the renormalized sequence x n := fi n / // n -i x t 



v z = J2^=o N{\z\ 2 ) z n I \Jx n \ e n . They also yield a resolution of the identity 

Jc N(\z\ 2 ) /i vr ^o Xnl 

and, by following the same procedure as above, a non-commutative version of the plane. 
For sequences which are perturbations of natural integers in a sense which will be made 

Kf{\z\ 2 ) 

of the standard case, with uncertainty relations modified as 



precise in the paper, the ratio „ ^(| 2 p) should be ~ 1. This represents a mild perturbation 



AQ AP = - (v z , {x N +i ~ x N )v z ) , 

where the diagonal operator xn is defined by xjy e n = x n e n . If the differences x n+ \ — x n 
are not too far from 1, then we deal with a kind of curved deformation of the square lattice 
of uncertainty cells corresponding to the standard case. 

Therefore, the aim of the paper is to explain how such sequences of numbers allow one to 
implement a specific coherent state quantization from the complex plane, different from the 
canonical quantization associated with the standard integers. Accordingly, consequences 
in terms of phase space localization are examined. 
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The rest of this article is organized as follows. In Section 2, we define what we call 
Delone sequences of nonnegative numbers, more precisely Delone a-perturbations of non- 
negative integers, and we give a list of examples of such sequences. We revisit in Section 
3 the quantization of the complex plane by using the ordinary (or Gaussian) coherent 
states. We then attempt in Section 4 to follow the same quantization procedure but this 
time based on an a-perturbation of the natural numbers. The existence of a POVM is 
dependent on the associated Stieltjes moment problem according to whether the latter 
does or does not have an explicit solution. Even when a solution exists, it is interesting 
to adopt an alternative strategy for constructing a POVM which, in a certain sense, is a 
perturbation of the standard coherent states POVM on the complex plane. This strat- 
egy is based on a "renormalization" of the original Delone a-perturbation of the natural 
numbers. We then formally proceed in Section 5 to the quantization of the complex plane 
arising from this renormalized sequence. We give in Section 6, some insights about the 
algebraic structure of the basic operators (ladder operators and the resulting successive 
commutators) emanating from this quantization. Next, we describe the localization as- 
pects of the quantization of the complex plane by examining the associated orthogonal 
polynomials generalizing the Hermite polynomials (Section 7) and the mean values (or 
"lower symbols") of operators obtained from this quantization (Section 8). Section 9 is 
devoted to a particularly illuminating example of a-perturbation, namely a situation in 
which generalized hyper geometric functions are involved and where the probabilistic as- 
pects of our approach are clearly brought out. In Section 10, we go back to the algebraic 
aspects of the quantization through a-perturbations of natural numbers by examining the 
associated generalized Heisenberg algebra, a concept developed by some of us. In Section 
11, we examine in a more mathematically oriented way the conditions under which it is 
possible to suitably Delone renormalize the a-perturbation of the natural numbers, and we 
present a new result on asymptotic Poisson and Gamma distributions. 

2 Delone sequences 

Suppose that we are given an infinite strictly increasing sequence of nonnegative real num- 
bers 



obeying the following two constraints : 

(dl) X is uniformly discrete on the positive real line ]R + : 3r > such that x n +i — x n > 
r for all n £ N, which means that there exists a minimal distance between two 
successive elements of the sequence, 

(d2) X is relatively dense on M + : 3 R > such that for all x £ K + 3 n S N such that 
\x — x n \ < R, which means that there exists a maximal distance, say L, between two 
successive elements of the sequence. 
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These conditions imply that lim^—xx, x n = oo. Actually, they mean that the sequence X 
forms a Delone set, and we will call X a (non-negative) Delone sequence in R + . We could 
choose the scale in such a way that x\ = 1 or L = 1, but we leave this as an option. 

The rationale behind this set of constraints is that the sequence X should thus appear 
as not very different from the set of natural numbers. It should be noticed that from this 
point of view, sequences like x n = n a , a 7^ 1, or x n = n logn are not Delone. Also, 
familiar deformations of integers, like g-deformations or (p, g)-deformations, 

N(p.9) = P p -i - q ' W« = ' (2) 

are not Delone. In this sense, these expressions could be viewed as singular deformations 
of the set N. 

To the sequence X corresponds the sequence of "factorials" x n \ = x\Xi---x n with 
xq\ = f 1, the "exponential" 



W) = E — > ( 3 ) 



n=0 Xn - 



and the sequence of "moment" integrals 



Xn^n = J o W) dt, (4) 

in which the appearance of the "corrective" factors ^ n 's is needed since there is no reason 
that the Stieltjes moment problem is solved for a generic pair {dt/M{t), x n \), as would be 
the case for (e~* dt,n\). 

The present study will more specifically concern sequences of numbers {x n } behaving 
asymptotically like natural numbers, x n ~ const, x n at large n. For convenience, we will 
put const. = 1 or equivalently we will consider the sequence x n /const.. We will impose 
more on this "natural" behavior, of course. We will require in our main example the 
following 

Definition 2.1 (Delone perturbation of N) Let N 3 n 1— > a(n) be a bounded function 
with values in the interval (—1, 1), ct(0) = 0, and such that its successive jumps a(n + 1) — 
a(n) have lower bound r — 1 with r £ (0, 1). Then the Delone sequence 

x n = n + a(n), n£N, (5) 

is called an a-perturbation of the natural numbers. 

Note that this r is a Delone lower bound for the sequence {x n } nS N in the sense of (dl). 

For instance, a(n) could be a constant shift, a(n) = e for n > 0. In order to fulfill the 
Delone condition (dl), we should then have —1 < e < 1 and < r < 1. Another simple 
case that will be considered in this paper is 

an + b 

0<nw a(n) = , (6) 

cn + d 
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with -1 < b/d < a/c < 1 or -1 < a/c < b/d < 1 and < r < 1. 
The perturbation function could also be periodic, like 



a(n) = e sinum , —1 < e < 1 , < r < 1 — 2\e sinw| . (7) 
The function a(n) could even be a random perturbation. 

Another interesting example is the set Zt" of non-negative beta-integers, i.e. all those 
positive real numbers which are polynomial in powers of an irrational real number (3 > 1, 
when they are expanded with respect to the base (3 using the usual greedy algorithm. 
When is endowed with specific properties (e.g. Pisot-Vijayaraghavan (PV) algebraic 
integers or more generally belonging to the so-called Parry family), these (3- integers form 
a quasiperiodic sequence with a finite number of possible adjacent differences x n _|_i — x n . 
The number of distinct distances between consecutive /3-integers is finite if and only if (3 is 
a Parry number. It is natural to ask in case of Parry numbers, whether the corresponding 
/3-integers behave similarly to the integers. Some answers concerning the additive, mul- 
tiplicative and diffractive properties of /3-integers already exist (see Gazeau(2006)| and 



references therein). In view of the "quasi" -periodic distribution of the /3-integers on the 
real line, it is pertinent to investigate the extent to which they differ from ordinary integers 
in an asymptotic sense, according to the nature of (3. 

Let us quote here precise results from |Balkova(2008)| that illustrate the similarity 
between sets N and Z^ = {b n \ n G N} for (3 being a Parry number by mentioning two 
properties: 

1. The limit cn := bim^oo — exists and belongs to the extension field Q(/3) on the 

n 

rational numbers. 

2. For j3 being moreover a PV number such that its minimal and its Parry polynomial 
coincide, the modulation a(n) = — — n is bounded. Parry polynomials are defined 
in |Balkova(2008)1 

Let us mention that both of the previous asymptotic characteristics are known for (3 
being a quadratic Parry unit. We underline that for quadratic numbers the notions of 
Parry and PV numbers coincide. The following proposition providing explicit formulae for 
/3-integers for a quadratic PV unit (3 comes from [Gazeau(2006)] . 



Proposition 2.1 If (3 is a simple quadratic PV unit, then 

b n = cpn + ^- + { 7-7-7; j ! n € n1 , where cp 



If (3 is a non- simple quadratic PV unit, then 



f M?}' neN }' where cp = 1 ~ W 
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where {x} denotes the fractional part of a nonnegative real number x. 



The simplest example is provided by the set Z+ of non-negative r-integers, where (5 is 
equal to the golden mean r = (l + \/5)/2. These r-integers form a quasiperiodic sequence 
with two possible adjacent differences x n+ \ — x n = 1 or 1/t. Their exact expression is 
given by 

= |b n = c T n- + |— | . n G N 1 , where c T = ~ 0.8541, 

where {x} designates the fractional part of a nonnegative real number x. Dividing by c T 
gives the "normalized" Delone sequence x n = b n /c T = n + a T (n) , with 



a T (n) 



n + 11 1 



1 + t 2 \ r 2 J r(l + r 2 
Since {x} S [0, 1), we deduce the bounds for the modulation: 

-0.1708 < cxrin) < - — — =r- « 0.2764. 



r(l + r 2 ) ' TV ' (1+r 2 

We will explain below how such Delone sequences of numbers allow one to implement 
a specific non-commutative (or "quantum") reading of the complex plane associated to 
Poisson-like probability distributions and to what extent this reading is "close" to the 
"canonical" reading provided by natural numbers. We will start with the sequence of 
natural numbers and then we will generalize the procedure to sequences X endowed with 
suitable properties. 



3 Natural numbers and the "fuzzy plane" 



Let Ti be a separable (complex) Hilbert space with scalar product (• , •) (antilinear on the 
left) and orthonormal basis eo, e\, . . . , e n , . . . , 



(fini Cn') — 0~nn' ■ 

To each complex number z € C corresponds the following vector in H : 



E 

n=0 



e 2 



(8) 



(9) 



Such a "continuous" set of vectors is called a family of standard coherent states in the 
quantum physics literature Klauder,Skagerstam(1985) Ali et al(2000)| . They enjoy the 
following properties. 
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(i) \\v z \\ = 1 (normalization). 

(ii) The map C 3 z \— ► v z is weakly continuous (continuity). 



(iii) The map N E n i— ► |(e n ,f^)| 2 = e z ' 2 — — is a Poisson probability distribution 

n! 

with average number of successes equal to \z\ 2 (discrete probabilistic content 
Ali,Gazeau,Heller (2008)1 ) • 

2 \ z \2n 

(iv) The map C 3 z i— >• I (e n , f 2 )| 2 = e z ' — — is a Gamma probability distribution (with 

respect to the square of the radial variable) with n as a shape parameter (continuous 
probabilistic content Ali,Gazeau,Heller (2008) | ) . 

(v) 

1=1 — P„ (10) 

where P z = (v z , -)v z is the orthogonal projector onto the vector v z and the integral 
should be understood in the weak sense (resolution of the unity in 7i). 



The proof of fjlQj) is straightforward and stems from the orthogonality of the Fourier expo- 
nentials and from the integral expression of the Gamma function which solves the moment 
problem for the factorial n\ : 



Ic 



— P z = V e n (e n ,,-)^L= / -e^Vf'=Ve„(e Bl .)=I. (11) 
Property (v) is crucial in the sense that it allows one to define: 

1. a normalized positive operator-valued measure (POVM) on the complex plane equipped 
with its Lebesgue measure and its a— algebra T of Borel sets : 

7T 

^Ah [ —P z €£(H) + , (12) 

where C(7i) + is the cone of positive bounded operators on 7i; 

2. a so-called Berezin-Klauder-Toeplitz quantization of the complex plane [Klauder(1963) 



Berezin(1975)| |Gazeau et al(2003)[ |Klauder(1995)] , which means that to a function 
omplex plane there corresponds the operator Af in TC defined by 

f~A s = —f(z,z)P s = Yl (Af) nn ,e n (e n ,,-), (13) 
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with matrix elements 

(^/)nn' = [ V /(*> , (14) 

provided that weak convergence holds. 

We emphasize here the point that we would like the reader to understand by a non- 
commutative reading of the complex plane |Madore(1995)| : a non-commutative algebra of 
operators in 7i takes the place of the commutative algebra of complex valued functions on 
the plane. Note that the map f|13j) is linear and, in view of (jlOp . the function f(z,z) = 1 
is mapped to the identity. These features are what we "minimally" expect from any 
quantization scheme. 

For the simplest functions f(z, z) = z and f(z, z) = z we obtain 

A z = a , ae n = \/ne n -i, aeo = 0, (lowering operator) (15) 
Az = a) , a) e n = \Jn + 1 e n +i (raising operator) . (16) 

These two basic operators obey the so-called canonical commutation rule : [a, a'] = I. 
The number operator N is defined by N = a) a and is such that its spectrum is exactly N 
with eigenvectors e n : Ne n = ne n . The linear span of the triple {a, a' ,1} equipped with 
the operator commutator [• , •] is the Weyl-Heisenberg Lie algebra, a key mathematical 
object for standard quantum mechanics. The fact that the complex plane has become 
non-commutative is apparent from the quantization of the real and imaginary parts of 

z = -j={q + ip) : 

A q d ^Q = -L(a + a^), A p d ^ P = ^L(a - a t) , [Q,P]=iI. (17) 
4 Generic Delone sequences and the "fuzzy plane" 

We now consider a generic Delone sequence as in ([T]). As a first remark, we easily assert 
the following result about the convergence of the associated exponential. 

Proposition 4.1 Let X = {£ n }neN be a non-negative sequence such that x n > for all 
n > 1 and linin^oo x n = oo. Then its associated exponential 



N{t) = Y,~r ( 18 ) 

has an infinite convergence radius. 



n=0 Xn ' 
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Proof. It is enough to apply the d'Alembert criterion for a power series of the form 

V°° a z n ■ 
Z^n=o Un/i ■ 

«n+l 



lim 

n— >oo 



lim = . 

n->oo x n+ i 



We then mimic the construction Q by considering the following family of vectors in 

H : 

00 1 z n 

" = Swwy73 e - (19) 

These vectors still enjoy some properties similar to the standard ones. 

(i) ||t> 2 || = 1 (normalization). 

(ii) The map C 3 z 1— > v z is weakly continuous (continuity). 

\z\ 2n 

(iii) The map N 6 n \— > \(e n ,v z )\ 2 = — - is a Poisson-like distribution in \z\ 2 , 

N{\z\ 2 )x n \ 

with average number of successes equal to \z\ 2 (discrete probabilistic content 



Ali,Gazeau,Heller(2008)[|Bardwell,Crow(1964)l ). 



\z\ 2n 



However, the map C 9 z 1— ► \(e n ,v z )\ = r — r is not a (Gamma-like) probability 

N(\z\ z )x n \ 

distribution (with respect to the square of the radial variable in the complex plane) with 
x^-j-i as a shape parameter, and this is a serious setback for the Berezin-Toeplitz quan- 
tization program. Indeed there is no reason why carrying out again the calculation (I13|) 
should yield the resolution of the unity : with P z = (v z ,-)v z , 

/^P.= £ «.<«,,.)-* I ^J^f* 

OO ^ 

= E— An)e n {e nr )=¥ (20) 

71=0 

where F is a diagonal operator determined by the sequence of integrals 

f + °° dt 

W e WY (21) 

These integrals form a sequence of Stieltjes moments for the measure 



We are now faced with the following alternatives. 
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The moment problem has a solution. More precisely, we first suppose that the 
Stieltjes moment problem has a solution for the sequence of factorials (x n \) n€N , i.e. there 
exists a probability distribution t \— > w(t) on [0, +00) |Ali,Gazeau, Heller (2008)] with infinite 
support such that 



rp I 
^71' 



t n w(t) dt . 

We know that a necessary and sufficient condition for this is that the two matrices 



/ 1 


Xl\ 


x 2 \ 


x n \\ 


/ Xl\ 


x 2 \ 


x 3 ! 




Xl! 


x 2 \ 


x 3 ! 


.. x n+1 \ 


x 2 \ 


x 3 \ 




• • x n+2 l 


x 2 \ 


x 3 l 


X4! 


• • x n+2 l 5 


X31 


X4! 


x 5 ! 


■ ■ x n+3 l 



\x n \ x n+x \ x n+2 \ 



x 2 J J 



\X n+ l\ x n+2 \ x n+3 \ . 



X2n+l } -J 



(22) 



(23) 



have strictly positive determinants for all n. 

Then, a natural approach is just to modify the measure in (|22[) by including the weight 
tt;(|z| 2 ) ftf(\z\ 2 ). We then obtain the resolution of the identity : 



/ ^(|z| 2 )AA(|z| 2 )P z = f e n (e n ,,.)^J== [ 



w(\z\ 2 )z n z T 



n,ra'=0 
00 



^ 1 r+00 

= J2~r t n w(t)dte n (e n ,-)=I. (24) 

n=0 Xn - J ° 

Of course, in many (if not most) cases the explicit form of the measure is not known. 
So, at this stage, the construction of CS may remain at a formal level and we will rather 
adopt the approach described below. 

The moment problem has no solution or has a solution with unsolved measure. 

Here we will suppose that the sequence of ratios 

(25) 

x n . 

has a finite limit as n — > 00, 

— ► [ioo < 00 ^ lim ^ n = 1 . (26) 

It is then natural to "renormalize" the non-negative Delone sequence ([TJ as follows. 

x n *== x n , n>l, x = , (27) 

Mn-1 
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and this entails 

S:J = — (xJ), neN. (28) 
Mo 



The assumption (|26p ensures that this renormalized sequence gets closer and closer to the 
original sequence, and the condition of the terms being strictly increasing is fulfilled beyond 
a certain point no- If some crossing happens for the first few elements of the sequence, 
we will reorder the sequence (j27|) so that it becomes strictly increasing; however, in order 
to avoid an inflation of notations, we shall maintain the same notation (possibly at the 
expense of some notational abuse). 

We then start again the construction in (|9|) by considering this time the renormalized 
family of vectors in TC : 

^ = £^^=-^=e n , (29) 

n=0 ^/M(\z\ 2 ) VXn - 

where Af(t) is the exponential associated to the new sequence, 

_ + °° +n 

N(t) = £ T7 ■ (30) 

n=0 

From (|26p it is clear that this series also has an infinite radius of convergence. Now, not 
only do the vectors (|29p enjoy properties similar to the vectors f)19[) . 

(*") 11^2 II = 1 (normalization). 

(it) The map C B z i— > v z is continuous (continuity). 

— \z\ 2n 

(Hi) The map Nenw \(e n , v z )\ 2 = — is a Poisson-like distribution with average 

N(\z\ 2 )x n \ 

number of occurrences equal to \z\ 2 (discrete probabilistic content), 

but they also enjoy what we need for carrying out a Berezin-Toeplitz like quantization, 
since : 

~ \z\ 2n 

(iv) The map C 3 z ^ \ (e n , v z )\ 2 = — is a (Gamma-like) probability distribution 

M(\z\ 2 )i n \ 

(with respect to the square of the radial variable) with x n +i as a shape parameter 
and with respect to the modified measure on the complex plane 

v def Kf(\z\ 2 ) d 2 z 

u(dz) = r 9 . (31) 

N(\z\ 2 ) /x vr 
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(v) The family of vectors (f29j) solves the unity : with P z = (v z , -)v z 
f ~ °° 1 f d 2 z 1 

v{dz)v z = £ e ^ e ™'>-)-7TrFT L 



^il(«)e„(e„,-) = I, (32) 



n,n'=0 
00 



ra=0 " 



the last equality resulting trivially from (|26|) and 

A statistical interpretation of the measure v(dz). As we already noted, the map 
N 3 n 1— > j-f is a Poisson-like distribution determined by the original sequence X = 
(x n ). On the other hand, the map 

NBiiH a (n) = — (33) 

/in 

defines the random variable cr(n). Hence, the ratio M {t) / hqN {t) is nothing else that the 
expectation value of a with respect to the Poisson-like distribution defined by X: 

Kf(t) ^ l l t n 



1 1 t n 



In consequence, the measure i/(cte) on the complex plane reads as 



K^) = (E^a)(|z| 2 ) — . (35) 



7T 



A comment regarding the "perturbed" Lebesgue measure v(dz) in (|31|) in the 
case where the moment problem has an explicit solution with resolving mea- 
sure w(t) dt : 



Even in this case, where an explicit solution to the moment problem (j22[) exists, it is 
interesting to follow the second route: keeping the perturbed Lebesgue measure (|31|) in the 
plane, carrying out the alternate quantization using it and comparing it with the quantiza- 
tion of the phase space C equipped with the non-symplectic measure w(\z\ 2 )J\f(\z\ 2 ) d 2 z/ir. 
We will examine in Section 18.21 this question in the simplest example of the constant shift 
a(n) = e, n > 1. We will also discuss in the conclusion the relevance of such different 
approaches to the "quantum processing" of the classical phase space, provided with the 
Lebesgue measure (for which classical states are just point or Dirac measures), or alterna- 
tively, provided with a different measure, giving a "mixed" nature to classical states. 
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Note also, that in the case which concerns us most in this paper, of a-perturbations of 
the natural numbers, since i — > 1, as n — > oo (see Eq. (f2o]) ). the discrete distri- 

bution, nH (AA(|z| 2 )) _1 |2;| 2ri /x n ! appears as a "perturbation" of the Poisson distribution. 
Similarly the continuous distribution |z| 2 \— ► (M\z\ 2 ))~ 1 \z\ 2n /x n \ can be looked upon as a 
"perturbed" Gamma-distribution. In this sense we also have here a "mild perturbation" 
on the original Bayesian duality, involving the Poisson and Gamma-distributions charac- 
terizing the case of the standard coherent states. 

5 Quantization of the classical phase space through renor- 
malized Delone sequences 

We know that Property (y) is crucial since it allows to define 

(vi) a normalized positive operator- valued measure (POVM) on the complex plane equipped 

d 2 z 

with its Lebesgue measure and its a— algebra T of Borel sets : 

7T 

^ A „/§!iQ^P l££(W)+ , (36) 
JaN(\z\ 2 ) /x vr 

where C{7i) + is the cone of positive bounded operators on 7i; 

(vii) the Berezin-Klauder-Toeplitz quantization of the complex plane : 

M{\z\ 2 ) d 2 z 

J{z,z)r z = 

n,n'=0 

where 



~ 1 f 1 d 2 z 

(A f ) nn , = / f(z,z)z n z r 

1 Z" 00 (in 



Jo fJ>oM{u) Jo 



\/x n \x n A Jo A t o-AA(n) 

(38) 

provided that weak convergence holds. 

For the simplest functions f(z,z) = z and f(z,z) = z we obtain lowering and raising 
operators respectively, 



A z = d , a e n = \/ x n e n -i , aeo = 0, (lowering operator) (39) 
A g = a* , a' e n = sj x n+ ie n+ i (raising operator) . (40) 
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These operators obey the non-canonical commutation rule : [5,5'] = xn+i — xn> where 
xyv is defined by xn = a) a and is such that its spectrum is exactly {x n , n G N} with 
eigenvectors xn e n = x n e n . The linear span of the triple {a,a' , Xn} is obviously not 
closed, in general, under commutation and the set of resulting commutators generically 
gives rise to an infinite dimensional Lie algebra. 

The next simplest function on the complex plane is the squared modulus of z, \z\ 2 = 
\{p 2 + <Z 2 )> i- e - one-half times the squared radius of the circle, in other words, the classical 
Hamiltonian of the harmonic oscillator when we consider the complex plane as the phase 
space of a particle moving on the line: 

A g g = a a 1 = x N+1 . (41) 

Therefore, the spectrum of the quantized version of \z\ 2 is the renormalized sequence (|27p 
up to a shift by 1. We will come back to this important fact in the conclusion. 



6 Localization properties and associated orthogonal polyno- 
mials 



Consider the operators: 



Q = -L[o + 5t], p = _L[a-6+]. (42) 



The operator Q acts on the basis vectors e n in the manner, 



n = V IT e ™~ 1 + V ~^2~ €n+1 ' ^ 

oo ^ 

Since for a Delone sequence the sum > diverges, the operator Q is essentially self- 



n=0 



adjoint [Borzov(2001)| and hence has a unique self-adjoint extension, which we again denote 
by Q. Let E\, A G 1, be the spectral family of Q, so that, 



Q= / A dE\ 



oo 



Thus there is a measure w(d\) on 1R such that on the Hilbert space L 2 (M, w), the action 
of Q is just a multiplication by A and the basis vectors e n can be represented by functions 
p n (A). Consequently, on this space, the relation (pl3|) assumes the form 



Apn(A) = c n p n _i(A) + c n+ ip n+ i{\) , c n = \ — , (44) 
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which is a two-term recursion relation, familiar from the theory of orthogonal polynomi- 
als. It follows that w{d\) = d(eo, E\eo), and the p n may be realized as the polynomials 
obtained by orthonormalizing the sequence of monomials 1, A, A 2 , A 3 , . . . , with respect to 
this measure (using a Gram-Schmidt procedure). Furthermore, for any ro-measurable set 
Aci, 



(e n ,E(A)e r 



PnWPmW tu(d\) 



and 



(Pn,Pm)i2(K,w) = / Pn(A)p m (A) w(d\) 



(45) 



(46) 



The polynomials p n are not monic polynomials, i.e., that the coefficient of A n in p n is 
not one. However, the renormalized polynomials 



are seen to satisfy the recursion relation 

q n +i(X) = A g n (A) - c 2 q n -i{X) , 



(47) 



(48) 



from which it is clear that these polynomials are indeed monic. Also, writing £ 2 for the 
vector N(\z\ 2 )z v z , when expressed as a function in L 2 (]R, ro), we see that the function 



G(z, A) := C(A) = E ^=f = E ^TT *»( A ) - 



(49) 



n=0 v ~" n=0 " 

is the generating function for the polynomials q n (or p n ), in the sense that 

1 x n \ d n G(z,X) 



9n(A) 



22 n! 



Si 



(50) 



2 = 



There is a simple way to compute the monic polynomials. To see this, note first that 
in virtue of (|43p and (|44j) . the operator Q is represented in the e n basis as the infinite 
tri-diagonal matrix, 

/0 ci ...\ 

ci c 2 
c 2 c 3 
c 3 c 4 
c 4 



Q 



(51) 
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Let Q n be the truncated matrix consisting of the first n rows and columns of Q and I n the 
n x n identity matrix. Then, 





( X 


-ci 




















\ 




-ci 


A 


-C2 
























-c 2 


A 


-C3 
























-C3 


A 


-C4 













AI n — Q n — 











-c 4 


A 






























A 


C n -2 
























• • -C„-2 


A 




-l 




^ o 

















Cn-l 


A 


/ 


now follows that q n is 


just the characteristic polynomial of Q 


ra ■ 












Qn 


(A) = 


det[AI n 


-Qn] ■ 









(52) 



(53) 

Indeed, expanding the determinant with respect to the last row, starting at the lower right 
corner, we easily get 



det[AI n - Q n ] = A det[AI 



n-l 



Q 



ra-U 



-"n—\ 



det[AI 



n-2 



Qn- 



(54) 



which is precisely the recursion relation (|48p . Consequently the roots of the polynomial q n 
(or p n ) are the eigenvalues of Q n . 



7 Localization properties of the quantum phase space viewed 
through lower symbols 

Given a function / on the complex plane, the resulting operator Af, if it exists, at least in 
a weak sense, acts on the Hilbert space 7i with orthonormal basis e n : the integral 

(V,!/V)= / f(z,z)\ty,v x )\ 2 u(dz) (55) 
Jc 

should be finite for all ip in some dense subset of 7i. It should be noted that if ip is 
normalized then (|55|) represents the mean value of the function / with respect to the 
V'-dependent probability distribution z \—> \(tfi,v z )\ 2 on the complex plane. 

In order to be more rigorous on this important point, let us adopt the following accep- 
tance criteria for a function to belong to the class of quantizable classical observables. 

Definition 7.1 A function C9zi-> f(z,z) G C is a CS quantizable classical observable 
via the map f i— > Aj defined by HFfy if the map C 3 z = |(g + ip) = (q,p) i— » (v z , AfV z ) is 
a smooth (i.e. £ C°°) function with respect to the (q,p) coordinates of the complex plane. 
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The function / is the upper |Lieb(1994) or contravariant |Berezin(1975) symbol of the 
operator Af, and (v z ,AfV z ) is the lower |Lieb(1994)| or covariant |Berezin(1975)] symbol 
of the operator Af. 

In Chakraborty et al(2008)] such a definition is extended to a class of distributions 
including tempered distributions. 

Hence, localization properties in the complex plane, from the point of view of the 
sequence {cc n }neN and its renormalized version {x n } n gN, should be examined through the 
shape (as functions of z) of the respective lower symbols 

Q(z) d ^ (v z , Qv z ) , P(z) ^ (v z , Pv z ) , (56) 

and the noncommutativity of the reading of the complex plane should be encoded in the 
behavior of the lower symbol (v z , [Q,P]v z ) of the commutator [Q,P]. The study, within 
the above framework, of the product of dispersions 

(a^Q) (a, Cz P) = i (v z , (x N +i ~ x N )v z ) = ^ (v z ,x' N v z ) (57) 

expressed in states v z should thus be relevant. Note that in the case of an a-perturbation of 
N, the r.h.s of this equality becomes \ (1 + (v z ,a'(N) v z ), where a n = x n — n. Introducing 
supremum and infimum i& for the set of "renormalized" fluctuations {a n+ i—a n , n G N}, 
we get the following bounds on the product of dispersions : 

X - (1 + i a ) < (A^Q) (A^P) < i (1 + s&) . (58) 



8 Delone sequences and generalized hypergeometric func- 
tions 

8.1 A statistical perturbation of the Bernoulli process 

It is well known that the Poisson distribution N 3 k i— ► e~ l t k /k\ with parameter t is the 
large n limit of the Bernouilli process of a sequence of n trials with two possible outcomes, 
one ("win", say) with probability p = t/n, the other one ("loss") with probability q = 1—p. 
The probability to get k wins after n trials is given by the binomial distribution : 

P k n) =^P k (l- P ) n - k , (59) 

and its Poisson limit is readily obtained. Now, imagine a sequence of n trials for which 
the probability to obtain k wins is given by the following "perturbation" of the binomial 
distribution : 

Pi n) = ■QnMp) , J>l n) (p) = l, (60) 

fe=0 
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where {x n , n G N} is an a-Delone perturbation of the natural numbers. Here, the sequence 
of unknown functions q n -k(p) should be such that, in the limit n — > oo with p = t/x n , we 
get our deformation of the Poisson law: 

i 1 t k 

i ' k i \ i i 

X n - k \x k l N(t) X k \ 

In order to have this limit, we first observe withp = t/x n that for an a-Delone perturbation 

2^n' h r \ X n X n —\ X n — k +\ t ( t \ t ( t \ fr>n\ 
P Qn;k{P) = T Qn;k ~> f ™ <?n;fc , (62) 



An) x n- k / \ IV foi\ 



X n _fc!xfc! ' X n X n X n X k \ ' \X n J n^oo x k \n^oo ' y X . 

and so we should have 

lim q n k ( — ) = -ttt-t • (63) 

n-»oo \X n J A/(t) 

8.2 The generalized hypergeometric function as an approximation of the 
exponential 

In view of Eq. [63l we would like to find an example of a perturbation which yields for 
the series M(t) something close to the ordinary exponential. A good candidate for this is 
precisely the confluent hypergeometric or Kummer function : 

ir/A „, , ,x q( a + l)...(q + n-l) t n 

N ft) := iiq(a a + e;t) = > r- -r ; — — 

7 ^ (a + e)(a + e+ 1) • • • (a + e + n - 1) ra! 



n=0 



t" 



n=0 Xn ' 

where 

x n = n + a(ra) = n H — p , (65) 

1 + — - 

n 

and where e is a small real parameter. We thus get an a-Delone perturbation of natural 
numbers which is of the type of the projective transformation ©. Moreover, we readily 
go smoothly back to the exponential as e — » 0. 

Now, from the integral representation of the Kummer function, 

T(a + e) f 1 .„ . Wl ± 



iJi(o; + e;t) = _. v ' ' / e tu u 8 " 1 (1 - u) 6 " 1 du , (66) 

r(a)r(e) y 

and from the usual polynomial approximation of the exponential, we derive the following 
binomial like approximation, 

iFi(a; a + e-t)= lim [ (l-—u) " * u 0-1 (1 - u)^ 1 du 



r(a)r(e) 7 V ^ 
lim 2 Fi(x n _ k ,a;a + e;t) . (67) 
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Hence, we get the explicit form for the functions q n ;k(p)' 



Qn;k(p) = — —FT? " : 7\ ( 68 ) 

W n 2 Fi{x n - k ,a;a + e;t) 



where the normalization factor — is defined by the condition 



£4%) = 1- (69) 

fc=0 



The simplest particular case of the x n given by (|65p . which gives a perturbation of the 
Bernoulli process, is for a = 1, i.e. x n = n + e, n > 1, xq = 0. In this case, we are able to 
solve the moment problem ()22[) with 



W (t) = r . (70) 

w r(e + i) y ' 

Now, following the second way by studying the sequence of moment integrals given by 

1 f°° t n dt 

" n = ^.L wr (71) 

with 

r(n + e + i) 
r(e + i) ' 1 j 

is straightforwardly calculated to be 

M(t) = iFi(l; 1 + e; t) = e i t" e [r(l + e) - e r(e, t)\ , (73) 
where r(e, t) is the incomplete Gamma function T(e,t) = J 4 °° x e_1 exp(— x)dx. Then 

As pointed out in the comment ending Section HJ we have here an elementary example 
allowing us to compare two types of non-commutaive reading of the complex plane. The 
first one stems from the solution of the Stieltjes moment problem and leads to the resolution 
of the identity with the measure 

,2> ,2><i 2 2 (, r(e,|2| 2 )\ Sz 



^T-l'-iRjT' (75) 

Note that in the limit z — > this measure reduces to (1 — e) ^f-. The second one is the 
measure associated with the renormalized sequence: 

u(dz) = ^ ^ = e~V ^ ^ (76) 

ndz) AT{\z\*)w r(i + e )- e r( e ,t) /.ovr' {fb) 
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where the AA(|z| 2 ) is numerically derived from the renormalized sequence (x n ). In Figure [T] 
we compare the two measures through the behavior of the two functions i i — >■ [ 1 



and t 



, r(M) 

noJUj^j ^ or ^ — 0- It is important to notice that the measure obtained by the 
solution of the moment problem implies a radical change in the statistical content of the 
complex plane as a phase space, since at t = the weight w(t)M(t) vanishes for e > 
and worse, it diverges for e < 0; in both cases it tends to 1 as t increases more or less 
rapidly depending on the value of e. On the other hand, the measure v issued from the 
renormalized sequence has a more reasonable behavior, whatever the sign of e, since it 
starts at l//xo at t = 0, remains strictly positive, and tends to l//Jo at large t, more or less 
rapidly depending on |e|. 




in 



n 
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Figure 1: The two measures for the sequence x n = n + e, n > 1, and xo = are compared 
through the behavior of the two functions t \— * (l — e ) (continuous line) and t i— ► 

JaJAt) nne ) f° r t > 0. The value of e has been chosen e = 0.5. 

A second example which is interesting to study is the particular case of ([6]) for which 
a = d = and b/c = — e, where e is a small real parameter. In that case, 



x n = n(l j ) ' n>l, x Q = 0, 



71- 



X 



and 



r(l-^ + n)r(l + y^ + n) 

r(i-v^)r(i + V£)n! 
Af(i) = 2^2(1,1;! - \/£,l + 



(77) 
(78) 
(79) 
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where 2-^2(^1) Q>2\ b\, 62; t) is a generalized hypergeometric function. 

As a matter of fact, the hypergeometric functions i-Fi(l;&;t) and 2-^2(1; 1; b±, 62; t) are 
particular cases of a more general class of functions which also yield perturbations of the 
exponential function, that is, q F q (li, I2, ■ b±, b 2 , ■■■bq; t), that are connected to Delone 
sequences of the form x n = n + a^ q \n), with 

a («) (n) = n? = i("-M-" g . (80) 

V > n q-l V > 



For these sequences, 



Yyi ^ r(l-6j+n) 



Xn: " pi=I ' (81) 
which gives for AA(i) exactly the generalized hypergeometric function, that is, 

Af(i) = q F q { M^^; 1 - 61, 1 - 6 2 , ...1 - b q ; t) . (82) 

q— times 

When q = 1 and 61 = —a, we have 

x n = n + a , n > 1 , xo = 0, (83) 

and 

A/"(i) = iFi(l;l+a;t) ; (84) 

for q = 2, 

(ra-6i)(n-6 2 ) 

x n = , n>l,x = 0, (85) 

n 

and 

AA(t)= 2^2(1,1; l-&i, I-62; t) . (86) 

The example ([77]) solved above is a g = 2 case, with a± = y/e and 02 = — \/e. 
There are sequences still more general than (|80p that result in perturbations of the 
exponential M(t), that is, 



Uli( n-h)-Ul=i (n-ai 
U q i=l( n - a i 



a W (n) = n llt=lv " I" , llt=l v " . (87) 



n ^ =1 ; ; , n>l,x = 0, (88) 



The Delone sequences for appropriate values of the constants {aj} and {bi} are then 

'Ui=i( n - a i) 

and our "perturbed exponential" becomes 

Af{t)= q Fq{l - ai ,l - a 2 , - a q ; 1 - h,l - b 2 , - b q ;t) . (89) 



The Delone sequence ([65]) is a particular case of ([88]) with q = 1, a± = 1 — a and 
61 = 1 — a — e. The constants {a^} and have to be chosen in such a way that the values 
of a^ q \n) satisfy the criteria required by the Delone sequences, for any value of n. 
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8.3 Asymptotic estimates 

Let us now analyze the asymptotic behavior of some sequences presented in the last para- 
grah for very large values of n. 

As a first example, let us take the Delone sequence (|65[) for the special value a = 1, i.e. 
the constant shift for n > 1. In (|74p . note that the numerator is the Gamma distribution 
which is centered around a value of t of the order of n. When we consider the limit for 
very large n, only large values of t contribute to the integral. But, for these values of t, 

e lM « i (90) 

r( e + i) << ' [ j 

and we have for \i n the approximate expression 

e Yin + 2e + 1) 

~ 1 + r(« + 1) r(n + \ + 2) fa. i - <■ * + * + * - i) . CD 

For large n, the hypergeometric function F(l, 1 — e; n + e + 2; —1) becomes a finite series; 
as n is much larger than e, we can keep only its first term. Therefore, we finally obtain 
that the asymptotic expression for fi n is 

— 1 + p7 r on - ( 92 ) 

which is pratically 1 for n = 10. In order to have a better idea of the perturbation in the 
spectrum x n introduced by fj, n , (see equation (JSTj)), we show in Figure [2] the behavior of 
x n — Xn for the first 15 values of n (e = 0.1). 

As a second example, let us consider the particular case of (JSJ) for which a = d = and 
b/c = —e, where e is a small real parameter. In this case, x n is given by the equation <\77\) . 
which leads to x n \ and M{t) shown by Equations (J7SJ) and (JT9"j) . fj, n can then be written as 

1 f°° — 

^"^Jio 2^2(1,1:1-^1 + ^;*) ' 

For large values of n, the main contribution for the integral of the previous equation is 
given by t ~ n and the dominant contribution of the hypergeometric for t >> 1 is 

2 F 2 (l,l;l-^,l + vM t> ~ i r[l- v^rfl + v^je* (l-^-fc^) , (94) 

leading to the following asymptotic behavior for \x n : 

Hn — ► 1 + 4 » ( 95 ) 
that presents the same behavior of x n /n in function of n as shown in Equation (|77p . 
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Figure 2: Delone sequence x n = n + e ,n > 1, and xq = 0. Behavior of x n — x n for the first 
15 values of n > (e = 0.1). 



A more complicated example of a Delone sequence with a decay in n similar to the first 
example of this subsection but presenting an oscillating character, is the sequence given by 
the following expression for a(n): 

esm(V2n) 

a(n) = , (96) 

n 

where e is a small parameter. In this !, M(t) and [i n can be calculated only 

numerically. We have found that the sequence fj, n oscillates, and the asymptotic behavior 
of its maxima is \fi n — 1| = e/n. In Figure [3] we show \x n for the first 300 values of n, taking 
e = 0.1 in (|96l). 



9 Sequences and triplets: a general algebraic setting 
9.1 Heisenberg-like triplet 

Suppose we are given a finite or infinite, sequence of nonnegative real numbers, bounded 
from below, {x n } n eN such that xq = 0. Let us view this sequence as the set of eigenvalues 
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Figure 3: Delone sequence x n = n + 0.1 , n > 1 , xq = 0. Behavior of // n , for the 

n 

first 300 values of n > 1. The large oscillations present near the origin have been truncated 
for convenience. 



for a positive self-adjoint operator denoted by x^. 

x N e n = x n e n . (97) 

The associated separable Hilbert space H is precisely the closure of the linear span of the 
orthonormal system {e n }, i.e. the latter is an orthonormal basis for 7i. We now define the 
operator labeling the states, or number operator, by N, 

Ne n = ne n , (98) 

providing also the rationale for the notation xn, or more generally of the notation f(N) 
for any diagonal real operator defined by f{N)e n = f{n)e n constructed from a sequence 
n — > f(n). The set M N of these sequences is a commutative algebra on K, and this naturally 
induces a structure of a real commutative algebra for the corresponding set of f(N). This 
algebra will be denoted by T>. We now define lowering and raising operators for the sequence 
Mn)}: 

o e n = \J x n e n —i , a e n = yj x n +\e n -\-\ . (99) 
From these definitions it is clear that 

x N = a)a. (100) 
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Now, for any diagonal operator f(N) in T>, the following commutation rules hold: 

of (TV) = /(AT + l) a = /(_i)(JV)a, f(N)a} = oV(-i)(JV), (101) 

[a,/(iV)] = /'(iV)a = a/^iV), [at,/(iV)] = -a^f'(N) = -f' {1) (N)a\ (102) 

where fr m -\(N) = f(N — m), m £ Z, is the m^ n shift of /, with the convention that f(n) = 
for n < 0, and /'(JV) = /(iV + 1) - f(N) is the finite difference derivative of f(N). Of 
course, the validity (or otherwise) of the above relations depend on domain considerations 
for the operators in question. In particular, we also have 

[a,at] = x' N . (103) 

The triplet (lowering a, raising a', commutator [a, a^] = x' N ) will be called a spectrum 
generating triplet of the sequence {x n }. In Gazeau,Lafortune,Winternitz(2008)] the Lie 



algebra generated by such a triplet has been investigated in a systematic way . 
9.2 Generalized Heisenberg triplet 

An algebraic structure which contains g-oscillators as a particular case was proposed in 



Curado,Rego-Monteiro(2001) and has been successfully used in some different physical 



problems Oliveira-Neto et al(2007) 


de Souza et al(2006) 


Bezerra,Curado,Rego-Monteiro(2002) 


Ribeiro-Silva,Curado,Rego-Monteiro(2008) 


Curado,Rego-Monteiro(2006) 


Curado et al(2008) . 



In this new algebra, called Generalized Heisenberg Algebra (GHA), the commutation re- 
lations among the operators a, and a diagonal operator Jo which is going to replace in 
the present context xn introduced in (jlOOp depend on a general functional h(x) and are 
given by: 

J a) = a^h(J ) (104) 

a,J = h(J )a (105) 

[a,<J] = h(J )-J , (106) 

with a = (at)t and Jo = Jq. This algebra tells us that the Jo eigenvalues ( Jo e n = 
t n e n ) are obtained by a one-step recurrence relation (t n = h(t n -i)), i.e, each eigenvalue 
depends on the previous one. Thus, the eigenvalue behavior can be studied by dynamical 
system techniques, simplifying the task of finding possible representations of the algebra 
|Curado,Rego-Monteiro(2001)| . In order to relate this algebra to a physical system, the Jo 
operator can be identified with the Hamiltonian of the system. 

When the functional h(x) is linear we recover Curado,Rego-Monteiro(2001)| the q- 



oscillator algebra; other types of functionals lead to other different algebraic structures 
|Curado et al(2001)] . The functional h(x) can be, for example, a polynomial and hence 
depend on some parameters (the polynomial coefficients). Let us consider the vector eo, 



2G 



with the lowest eigenvalue of Jo, Jo eo = toeo- For a general h(x), these operators act on 
the Fock space as 

J e m = /i (m) (t )e m , m = 0,1,2,- ■■ , (107) 
a f e m = N m e m+ i, (108) 
ae m = iV m _ie m _i, (109) 

where N% l _ l = (to) — to and /i^ m ^ (to) is the m-th iteration of to by means of the function 
h. 

It is interesting to observe that in some cases we can find a GHA associated to Delone 
perturbations of N. Indeed, let us consider those Delone sequences, eq. ©, that can be 
inverted. In those cases we can write 

n = x n + j(x n ) . (110) 

Note that the difference (x n — n) implies that ")(x n ) shares the same properties as a(n). 
Therefore, the inverse sequence is also a Delone sequence. 

In order to find the GHA for invertible Delone sequences, we write the (n + l)-st term 
of the sequence in terms of x n : 

x n+ i = x n + 1 + 5(x n ) , (111) 

where 

5(x n ) = 7(cc n ) + a(x n + 1 + j(x n )) . (112) 
The characteristic function is trivially obtained from eq. (jl 121) : 

h(x) = x + l + 5(x) , (113) 

and the corresponding GHA is then 

[Jo,A^] = A\l + 5(Jo)) (114) 

[a,a f ] = 1 + 6(J Q ). (115) 

From eqs. ([5]) and (jllOp it can be seen that a and 7 are perturbations of opposite signs. 
Thus 5 is the difference of two perturbation terms and hence itself a perturbation. Eqs. 
(|114p and (|115p show that the invertible Delone sequences associated to GHA are pertur- 
bations of the quantum harmonic oscillator. 

We take as an example of an invertible Delone sequence the case 

CLTl 

a(n) = -, n>0, d^0, (116) 

cn + a 

where the real constants a, c and d have to satisfy the following conditions: 



4 

< 1; c^0; — ? ad > r - 1 , (117) 



c 2 
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Figure 4: Delone sequence x n = n + a(ra) with a(n) = c °? d , for the case a = 0.1, c = 1, and 
d = 2.2. The behavior of a(n) (dashed line), 5(x n ) (continuous line) and "f(x n ) (dot-dashed 
line) is shown for the first 10 values of n. 



and r 6 (0, 1), according to Definition 2.1. 
For this sequence, 



J c 2 x\ + 2c(d - a)x n + (a + cf) 2 - cx n - a - d 

7(»n) = J" ( 8 ) 

and 

For the particular case x n = n + a{n) with a(n) = e ^" rf , see Figure H] where the 
functions a, 7 and 5 are explicitly exhibited. 

If we choose specifically a = e, |e| < 1, c = 1 and c£ = 0, then a(n) = e, 5 = and the 
GHA is given by h(x) = x + 1 which would characterize the harmonic oscillator algebra. 
However this sequence is out the scope of the present paper since xq = e 7^ 0. Imposing 
xq = besides x n = n + e for n > 1 takes us to the case already examined in the previous 
section and for which we do not have a GHA algebraic structure. 



10 A result on asymptotic Poisson and Gamma distributions 

Let us now present a result which explains the behavior of the sequence (fi n ) given a certain 
class of Delone sequences, namely the existence of fj,^ for o-Delone perturbations of N. 
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For an a-Delone perturbation x n = n + a(n) of N one can write 

1 1 «Tt/i, 



x n ; = n! 



««) «(«) 



Then the associated "exponential" reads as 



00 y.71 

^— ' n! 

n=0 

Its ratio to the ordinary exponential reads as the Poisson average of the random variable 

n i — > 

e ~ 4 = £ ^ e ~^) = ( E ^) (*) > 

n=0 

Thus, the ratio fJ> n /£,( n ) can be rewritten as the Gamma average of the random variable 
1 1 — ► e'/.A/"(t) : 



«n) f(")Jo V(t)x„! i„ n! jV(t) V u (E P f)(t) 

Let us put this on a more general setting. Let N 9 n 1— ► £(n) be a discrete function 
which is extendable to a function R 9 cc 1 — ► M with £(0) = 1 and lim^—xxj f (x) = 0. Its 
Poisson mean value with parameter t is given by 

(E P 0(t):=E^ e ~^)' ( 121 ) 

n=0 

whereas the Gamma mean value with parameter m of a random variable S is given by 

(E G H)(m)=/ -r e-*S(t)dt. (122) 
Let us examine the asymptotic behavior of the following combination of these two averages: 



t 



rn 



1 \ r°° - r^( m ) 

I m = £(m) ( E G — — (m) = / m! , re (it . (123) 



(EpO; 7o E - Q ^(n) 

We would like to give sufficient conditions for having linim^oo I m = 1. 
Proposition 10.1 Let us define the logarithm o/l/£ by <p : 

£(x) = e -*(«) , (124) 
and suppose that the function IR + 9 x 1— ► </>(x) /ias i/ie following properties. 
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with f (a) = and 



(%) is three-times derivable almost everywhere (a.e.) on ~R + , 

(ii) \(p'(x)\ < lfa.ej ai Zarye sc. ^Zso x|0"(sc)| <C l(a.e.) and x 2 \(p'"(x)\ < lfa.e.j at 
Zarge x. 

T/ien I m — > 1 . 

m— >oo 

Proof. 

Let us put 

oo n oo 

F t (x) := xln- + x - -Inx - 4>(x) (125) 
3/ 2 

where we have neglected the errors in the Stirling's formula, n n e _n (\/2vrn)/n! = l + 0{n~ l ) 
(e.g. e>(?i _1 ) = -0.012 for n = 7.) It then follows 

Fi(x) = In * - 1- - tf/(x) , if (x) = -1 + -L _ ^( X ) . 
x 2x x 2ar 

We assume that there exists a solution to F((x) = at x = x*(t). This implies 

x * e ^+<t>'{*') = t ( 12 6) 

We now impose two conditions. 

Requirement l^We demand that [t — > oo/ implies [x*(t) exists and that x*(t) — ► oo] 

In order for x*(t) to exist at large i, the function xe2^ +<?!> ^ needs to be unbounded. 
We face two cases: 

• (Case-a) If it diverges for finite x due to a singular behavior of <f>'(x), then the second 
condition (x*(t) — > oo) is violated. Therefore this case is excluded, i.e. <p'(x) should 
be locally bounded. 

• (Case-b) The asymptotic behavior of (fr(x) for a; > 1 is such that it allows the 
divergence of xe^ + ^ ^ in this limit, i.eH <fi f {x) can be negatively large but it does 
not dominate Inx. That is, either \4>'{x) > ] or [4>'(x) < and |0'(x)| <C Inx] for 
x > 1. 



1 This requirement comes from the application of the saddle point method. 
2 xe^ ~ (x+ §) for x > 1. 
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In summary 

<!>'{%) > — lnx + (positive increasing) (for x — > oo) 
The condition for </>(x) in the hypothesis is modest enough to satisfy the above condition. 
Requirement 2JI For x* 1, the second derivative at the peak \F"(x*)\ is negative, and 

\F?{x*)\-* «x* 

so that the hump of the exponential becomes asymptotically sharp. The above inequality is 
equivalent to (assuming that the negative value condition is satisfied )@ 

1 1 

~ + 4>"{x) > -5. 

X X L 

The second derivative 4>"(x) can be negative, but does not dominate or cancel the -. We 
formulate this condition as 

4>"{x) > 1- r](x), rf(x) S> — iy (for x — > 00) 

x x z 

The condition for (f)(x) in the hypothesis is modest enough to satisfy the above condition. 



Application of the Laplace method: 

With the above two requirements being satisfied, the evaluation of the sum is done by 
using 

1 

wf *\ \Fl'(x*)\ (x-x*) 2 

e Ft(x) ^ e Ft(x*) e 2 . 



The result is H 

A(t) 



Note that 1 < y/x*. 

y/l+x*<t>"(x*) 

Next we evaluate 



exp(i + x* <j>'(x*) + x* - ct){x*)) 
^l + x*<t)"{x*) 



DC 



with 



Gr(t) = mln h m — 0(m) — In A(t) 

m 



3 This requirement also comes from the application of the saddle point method. 

4 We can neglect — with respect to i 

5 We used x* \F{'(x*)\ ~ 1 + ^"(a;*) for a:* > 1. 
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~ m In h (m — x*) — [(j)(m) — <j){x*y\ 

m 

--- x*4>'{x*) + - ln(l + x*(j>"{x*)) 
The time derivative of the defining equation of x* (eq J126p gives 



e2x 



dx* 



dt 1 + x*(j)"(x*) 



(for t — > oo) 



Then 



G'{t) 



m dx* 
T It 



-(l + xW)) + 



+ 





>"'(x*) 


2(1 + x*</>" 


^*)) 


+ x *(j>'"(x*y 


x*4>"{x*)) 2 _ 





(127) 



(128) 



Note that the denominator in [] is such that (1 + x*(j)"(x*))~ 2 < x*. 
Final concession: 

Here we shall introduce further simplifying assumptions ( at the possible cost of narrow- 
ing the validity region of JH-" 

x*<t>"{x*) + x* 2 cf)'"{x*) -► (for t -> oo) 
1 



2x* 



+ 0'(x*)^O (fort^oo) 



These conditions are just fulfilled under the hypotheses of the proposition. 



Second application of the Laplace method: 

Then we have first x* ~ t asymptotically, and G'(t) ~ ™ 

G'{m) = (for t oo) 



1. Therefore, 



and G"(t) = -p gives 

The Laplace method then yields 
I ~ 1 

J- 777. 



V2 



irm Jo 



G"(m) 



1 

m 



V2 



irm 



e 2m dt = 1. 



If we wish to improve the validity of the theorem, we should put more precise conditions here. 
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Implications : 

As for £(m) = e~^ m \ the conditions 

4>'(x) -> 0, ^"(x) 0, x 2 (j)'"{x) -► (for t oo) 

are rather loose. 4>{x) = {—a) lnx with both a > and a < are permissible. That is £(m) 
can grow or decrease as a power law, £(m) ~ x a . So 4>{x) = ax q with < 1 and both a > 
and a < are admissible. That is, £(m) can even grow like e am<J with a > and g < 1. 

11 Application to Delone a-perturbations of N 

Let us consider an a-perturbation of N 

x n =n + a(n), n G N , (129) 

where N3iih a(n) is a bounded function with values in the interval (—1, 1), a(0) = 0, 
and such that its successive jumps a(n + 1) — a(n) have lower bound r — 1 with r £ (0, 1). 
One can write x n l = n\ -J-^ with 

J- = ft + 



k=l 

Then the logarithm of l/£ reads as : 



k=i v J 



(130) 



Let us extend the discrete domain of this function to a continuous one by replacing the 
sum by the integral 



Hence, 



</>(x) = J\(l + ^\dy. (131) 

<P'(x)=ln(l + ^-) « ^ (132) 
V ^ / at large x a; 

and so </>'(x) — ►x-xxj since x t— > a(x) is bounded. For the second derivative, we have 

„ xa'Qr) - a(g) a'(a;) 

a; 2 + xa(x) a t large x x 

Therefore, from Proposition 110.11 one can assert the following. 
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Proposition 11.1 Let x n = n + a(n) , n G N , an a-perturbation of N such that the map 
N 3 n a(n) extends to a function R + B x i— > a(x) £/iaf is a.e. two times derivable. 
Suppose that \a'(x)\ <C 1 and x\a"(x)\ <C 1 at large x. Then we have 

Im = / n dt ^ m ^oo 1 . (134) 
JO y^oo _E 

An example of such a sequence is given by the periodic function : 

o(x) = \{nx} + v , (135) 
with |A| < , r G (0, 1), |a'(x)| = |A/x| a.e., a"(x) = a.e., and 

-1< v < 1 - A if A>0, 
-A - 1 < i/ < 1 if A<0. 

Note that the sequences Z^/cg of rescaled non-negative /3-integers with /3 a quadratic PV 
unit, such they are defined in Proposition 12.11 are of this type and so fulfill the conditions 
of Proposition 111.11 



12 Discussion and outlook 

In this work we have explored the possibility of quantizing the complex plane through the 
use of sequences of numbers close enough to the natural numbers. We have discussed the 
probabilistic content of the procedure. Some algebraic aspects have also been considered. 
We have treated a specific example, the hypergeometric case, which leads to explicit esti- 
mates. A general mathematical result has also been presented which pertains to the present 
context as well as to mathematical statistics. Many aspects of our work deserve further 
investigation, on different levels, numerical, mathematical, and interpretational, specially 
those around the notion of quantum or fuzzy localization in the complex plane arising from 
such a "non-commutative" reading of the complex plane. 

Let us just discuss, from a more oriented physical point of view, the relations between 
the original sequence X = {x n , n G N}, its renormalized companion X = {x n , n G N} 

Kf(\z\ 2 ) d 2 z 

and the modified Lebesgue measure on the complex plane, v(dz) = r . Suppose 

M{\z\ z ) /i vr 

we observe through some experimental device the sequence of numbers X, for instance 
the quantum energy spectrum of a given system, that lies in the class of a-perturbation 
of the nonnegative integers. Due to this proximity, probably the considered system is 
classically described by a harmonic oscillator with Hamiltonian H = + q 2 ). Since 
the canonical quantization of the harmonic oscillator yields the integers as a spectrum, we 
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here attempted to set up a quantization framework which takes account of the modified 
spectrum. For this purpose, we need to solve a Stieltjes moment problem. We know that 
this is impossible for most of sequences X, for instance the sequences of rescaled beta- 
integers mentioned in this paper. Hence we are naturally led to deal with the renormalized 

sequence X for which we solve the moment problem with the measure -rr, \ — • This 

N{t) no 

measure endows the complex plane viewed as the phase space for the harmonic oscillator 
with a statistical content, in the sense that classical states are not anymore described by 
points zo, i. e. by Dirac distributions z \— > 5 Zo (z), but instead by a smooth distribution of 

the type z \— > Z °\ J . Then, by following the CS quantization stemming from the 

■A/ O - zon 

renormalized sequence, we consistently find the latter as the renormalized version of the 
observed spectrum along the equation 



KP 2 +9 2 ) 



A z2 = aa* = x N+ i 



This spectrum is not the observed one, but it might differ appreciably of the later in the 
first levels only, since x n = (fi n /fi n -i) x n . Actually only the differences x n+m — x n are 
experimentally significant and we should reasonably expect that they are practically the 
same as x^-j-^ x^i for excited levels, as it can be observed in the particular example shown 
in Figure [2j 



Acknowledgements The authors are indebted to B. Heller for fruitful comments and 
bibliographic suggestions. J.P.G. acknowledges the Centro Brasileiro de Pesquisas Ffsicas 
(CBPF) for financial support. 



References 

[Ali et al(2000)] Ali, S. T., Antoine, J. -P., & Gazeau, J. -P., Coherent states, wavelets and 
their generalizations, Graduate Texts in Contemporary Physics, Springer- Verlag (New 
York, 2000). 

[Ali,Gazeau,Heller(2008)] Ali, S. T., Heller, B., & Gazeau, J.-P., Coherent States and 
Bayesian Duality, 2008, J. Phys. A: Math. Theor., 41, 365302 1-22. 

[Balkova(2008)] Balkova, L., Gazeau, J. P., & Pelantova, E., Asymptotic behavior of beta- 
integers, 2008, Lett. Math. Phys., 84, 179-198. 

[Bardwell,Crow(1964)] Bardwell, G. E., & Crow, E. L., A Two-Parameter of Hyper-Poisson 
Distributions, 1964, J. Am. Stat. Association, 59, 133-141. 

[Berezin(1975)] Berezin, F. A., General concept of quantization, 1975, Comm. Math. Phys., 
40, 153-174. 



35 



[Bezerra,Curado,Rego-Monteiro(2002)] Bezerra , V. B., Curado, E. M. F., & Rego- 
Monteiro, M. A. , Perturbative computation in a generalized quantum field theory, 
Phys. Rev. D, 66, 085013 1-9. 

[Borzov(2001)] Borzov, V.V., Orthogonal polynomials and generalized oscillator algebras, 
Integral Transforms Spec. Funct., 12, 115-138. 

[Chakraborty et al(2008)] Chakraborty, B., Gazeau, J. P., <fe Youssef, A., Standard coherent 
state quantization of functions and distributions, submitted. 

[Curado,Rego-Monteiro(2001)] Curado, E. M. F., & Rego-Monteiro, M. A., Multi- 
parametric deformed Heisenberg algebra: a route to complexity, J. Phys. A: Math. 
Gen., 34, 3253-3264. 

[Curado et al(2001)] Curado, E. M. F., Rego-Monteiro, M. A., & Nazareno, Heisenberg- 
type structure of one-dimensional quantum Hamiltonian, H. N., Phys. Rev. A, 64, 
012105 1-5. 

[Curado,Rego-Monteiro(2006)] Curado, E. M. F., Rego-Monteiro, M. A., Rodrigues L. M. 
C. S., & Hassouni, Y., Coherent states for a degenerate system: the Hydrogen atom, 
Physica A, 371, 16-19. 

[Curado et al(2008)] Curado, E. M. F., Hassouni, Y., Rego-Monteiro, M. A., & Rodrigues, 
L. M. C. S., Generalized Heisenberg algebra and algebraic methods: the example of 
an infinite square-well potential, Phys. Lett. A, 372, 3350-3355. 

[de Souza et al(2006)] de Souza, J., Oliveira-Neto, N. M., &; Ribeiro-Silva, C. I., A method 
based on a nonlinear generalized Heisenberg algebra to study the molecular vibrational 
spectrum, Eur. Phys. J. D, 40, 205-210. 

[Erdelyi(1956)] A. Erdelyi, Asymptotic expansions, Dover, New York (1956). 

[Gazeau(2006)] Gazeau, J. P., Verger-Gaugry, J. L.: Diffraction spectra of weighted Delone 
sets on beta-lattices with beta a quadratic unitary Pisot number, Annales de l'lnstitut 
Fourier, 56, 2437-2461. 

[Gazeau et al(2003)] Gazeau, J-P. , Garidi, T., Huguet, E., Lachieze-Rey, M., & Renaud, 
J., Examples of Berezin-Toeplitz quantization: Finite sets and Unit Interval, 2003, 
Symmetry in Physics. In memory of Robert T. Sharp, Ed. Winternitz, P., (Montreal: 
CRM Proceedings and Lecture Notes 2004) . 

[Gazeau, Klauder(1999)] Gazeau, J-P. , Klauder, J., Coherent States for Systems with Dis- 
crete and Continuous Spectrum, 1999 J. Phys. A: Math. Gen., 32, 123-132. 

[Gazeau, Lafortune,Winternitz(2008)] Gazeau, J-P. , Lafortune, S., & Winternitz, P., Al- 
gebras generated by spectrum generating triplets, in preparation. 



36 



[Klauder(1963)] Klauder, J. R., Continuous-Representation Theory I. Postulates of 
continuous-representation theory, 1963, J. Math. Phys., 4, 1055-1058; Continuous- 
Representation Theory II. Generalized relation between quantum and classical dy- 
namics, 1963, J. Math. Phys., 4, 1058-1073. 

[Klauder, Skagerstam(1985)] Klauder, J. R., & Skagerstam, B. S., Coherent States - Appli- 
cations in Physics and Mathematical Physics, World Scientific, (Singapore, 1985). 

[Klauder(1995)] Klauder, J. R., Quantization without Quantization, 1995, Ann. of Phys., 
237, 147-160. 

[Lieb(1994)] See Lieb, E., in Feng, D. H., Klauder, J. R., & M. Strayer (eds.) : Coherent 
States: Past, Present and Future, World Scientific, (Singapore, 1994). 

[Madore(1995)] Madore, J., An Introduction to Noncommutative Differential Geometry and 
its Physical Applications, Cambridge University Press (Cambridge, 1995). 

[Magnus et al(1966)] Magnus, W., Oberhettinger, F., & Soni, R. P., Formulas and The- 
orems for the Special Functions of Mathematical Physics, 3rd ed., Springer- Verlag, 
(Berlin, Heidelberg and New York, 1966). 

[Oliveira-Neto et al(2007)] Oliveira-Neto, N. M., Curado, E. M. F., Nobre, F. D., & M. A. 
Rego-Monteiro, A simple model to describe low-temperature behaviour of some atoms 
and molecules: an application to the Hydrogen atom, 2007, J. Phys. B, 40, 1975-1989. 

[Ribeiro-Silva,Curado,Rego-Monteiro(2008)] Ribeiro-Silva, C. I., Curado, E. M. F. , & 
Rego-Monteiro, M. A., Deformed scalar quantum electrodynamics as a phenomeno- 
logical model for composite scalar particles, 2008, J. Phys. A: Math. Theor., 41, 145404 
1-18. 



37 



